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ABSTRACT 

Following previous suggestions of other researchers, this paper discusses the 
prospects for astrometric observation of MACHO gravitational microlensing events. 
We derive the expected astrometric observables for a simple microlensing event with 
either a dark or self-luminous lens, and demonstrate that accurate astrometry can 
determine the lens mass, distance, and proper motion in a very general fashion. 
In particular we argue that in limited circumstances ground-based, narrow-angle 
differential astrometric techniques are sufficient to measure the lens mass directly, 
and other lens properties (distance, transverse motion) by applying an independent 
model for the source distance and motion. We investigate the sensitivity of differential 
astrometry in determining lens parameters by Monte Carlo methods, and derive a 
quasi-empirical relationship between astrometric accuracy and mass uncertainty. 



Subject headings: gravitational microlensing, astrometry, dark matter 



1. Introduction 



In 1986 Paczyhski ( Paczyhski 1986b| ) suggested that photometric observations of gravitational 



microlensing might be used to indirectly study the population of massive compact objects in the 
galaxy, and in particular MAssive Compact Halo Objects (MACHOs) that might be a significant 
component of the dark matter thought to exist in the galaxy by dynamical considerations. 
Paczyhski's 1986 paper, and the observational proposals it fostered were met with some skepticism. 
However, the past several years have seen Paczyhski's suggestion spectacularly confirmed - at 
the time of this writing four separate groups have reported significant numbers of candidate 
gravitational microlensing events from photometric observations of LMC, SMC, and galactic bulge 
sources. The large majority of the light curves for these candidate microlensing events match 
theoretical expectations for single lens objects, and all collaborations report a significant excess 
of microlensing event candidates above the number expected from known stellar populations. 
In particular, from their first two years of data the MACHO collaboration reports eight events 
toward the LMC where only one is expected from known stellar populations, and estimates that 
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roughly half of the expected dark matter in the galactic halo is in the form of dark stellar mass 
objects (Alcock et al 1996). The difficulty in interpreting the MACHO collaboration events is 
that they are observed photometrically, which does not uniquely determine the mass of the lens 
- instead the MACHO collaboration bases their conclusions on interpreting their event sample 
observables (namely event duration) in the context of a halo model (Alcock et al 1996). However, 
other interpretations are possible (see |Sahu 1991 , Zhao 1997 , and |Gates et al 1997 for recent work 
arguing against a halo interpretation). 

Clearly it is desirable to measure MACHO physical properties in a model-free context. 
This objective has led a number of authors to propose the astrometric observation of MACHO 
gravitational microlensing events ( |H0g et al 1995 , Miyamoto fc Yoshi 1995 , Walker 1995| ), a 
specialized application of an earlier suggestion by Hosokawa et al ( posokawa et al 1993[ ). In 
particular, Miyamoto and Yoshi proposed the separate astrometric observation of both lensing 
images in MACHO microlensing events (a small misuse of the term microlensing - see [Paczyriskl 



1986a), and developed the theory of such astrometry. As we shall argue below, we find this 



suggestion implausible because of the small separation of the images. Instead, herein we consider 
astrometry of the lensed center-of- light, primarily for dark lenses. We find, as did Miyamoto 
and Yoshi, that high-precision astrometric observation of such microlensing events allows the 
estimation of the lens parameters (mass, distance, proper motion) appealing only to the properties 
of lensing. Moreover, we find that in a limited set of circumstances, the problem of determining 
a subset the lens parameters (mass, relative parallax, relative proper motion) is amenable to 
narrow-angle differential astrometric techniques very similar to those proposed and employed in 
gravitational companion search programs ( [Shao &: Colavita 1992 , Lestrade et al 1994, Benedict et 



al 1995, Gatewood 1996). In particular we argue that if a suitable astrometric reference frame can 
be established, the lens mass can be directly measured by ground-based differential astrometric 
techniques independent of additional assumptions, and the lens distance and transverse velocity 
can be estimated by appealing to an independent model of source distance and proper motion (see 
similar remarks by |Walker 1995| ). In such circumstances, many of the current issues regarding the 
nature of the lensing objects can be resolved. 

In this paper we assess the ability of astrometry to probe the physical parameters of 
microlensing events in which the lens is dark, with a particular emphasis on MACHO microlensing 
events. In §^ we introduce the theory to analyze a microlensing encounter as observed by a 
(near) terrestrial instrument, in terms particularly oriented toward narrow-angle differential 
astrometry. In §|| we address astrometric sensitivity to microlensing parameters through Monte 
Carlo techniques. Finally, in §|I] we place our results in the context of envisioned astrometric 
instrumentation, discuss the near term prospects for such an astrometric program, and mention 
future extensions to this work. 



2. Microlensing Encounter Description 
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2.1. Instantaneous Theory 



As quantitatively described by Refsdal (Refsdal 1964) and recently reviewed by Paczyrisk 



(Paczynski 1996a), a gravitational lensing event is photometrically observable when an intervening 
compact massive object passes close to the line-of-sight to a background source. This phenomenon 
follows from the curvature of space near a massive object. General Relativity predicts that an 
object of mass m deflects a light ray by an angle a: 

AGm 

a 



d 2 b 



where b is the "impact parameter" (transverse separation at the point of closest approach) of the 
light ray relative to the mass position. In this way the mass acts as an optical lens. 

The instantaneous geometry of a background source lensed by an intervening point mass 



(the lens) is depicted in Figures 1 and 2 of Paczynski 1996a. When the lens is sufficiently close 



to the nominal line-of-sight to the background source, an observer (equipped with a telescope 
of arbitrarily high angular resolution) sees a background disk-like source as two arcs at distinct 
positions corresponding to two solutions of a quadratic equation in the bend angles a, - 
on opposite sides of the lens. In a plane normal to the unperturbed line-of-sight to the source 
containing the lens (herein referred to as the transverse lens plane), the apparent positions of the 
two arcs are at impact parameters b± and 62 relative to the lens, functions of the impact parameter 
of the unperturbed line-of-sight to the source b s relative to the lens: 

b 1>2 = \ (±b s + ^ + 4i?|) = ^ (±u + v/^TI) (1) 

where all the impact parameters are positive semi-definite quantities, the subscript 1 refers to the 
positive quadratic root - the brighter of the two images and the closest to the source (in contrast 



to |Paczyhski 1996a| ), Re is the so-called Einstein radius: 



_ AGm D sl D lo UGm 
R E = \l— 2 ^ = ]l-^D so x(l-x) (2) 



and we have introduced a dimensionless impact parameter u = b s jRjr;. Dsh -^loi and D so are 
the source-lens, lens-observer, and source-observer distances, respectively, and x is the fractional 
separation x = D[ /D so . The two images appear separated in the transverse lens plane by: 



Ab = h + b 2 = yJV* + 4Jt% = R E Vu 2 + 4 > 2R E (3) 

The intensities of the two images can be computed from the fact that any lensing conserves 
surface brightness. For imperfect lens alignment and constant surface brightness of the source the 
relative brightness between the arcs and source is given by the ratio of perturbed and unperturbed 



image surface areas in the transverse lens plane of Paczynski 1996a Figure 2. The relative 
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intensities (units of the unperturbed source intensity) of the two images of a quasi-point source 



are easily shown to be (Paczyhski 1996a): 

u 2 + 2 1 



11,2 



2uVu 2 + 4 2 



±« (4) 



If we consider lensing of LMC, SMC, or bulge stars by intervening star-like compact masses, 
then what we observe is limited by the finite resolution of our instrumentation. In particular, 
for m ~ Mq and Di ~ D so /2 [x ~ 1/2) we find a value of Re (hence Ah) that is a few AU. 
For a lens distance on the order of kiloparsecs, that makes the images separated by an angle on 
the order of a milliarcsecond (mas - 10 -3 arcseconds). This image separation is well below the 
angular resolution of available instrumentation, so an observer sees the lensed images of the source 
as unresolved and of total amplitude ( [Paczyhski 1996a ) : 



ii 2 4- 2 

A = A 1+ A 2 = ^— >1 (5) 
uyu + 4 



The fact that the source images are unresolved qualifies this event as microlensing (Paczyhski 



1986a, Paczyhski 1996a| ). It is important to note that A is significantly greater than one only for 



u 



less than one (Figure [TJ) , so the unperturbed source line-of-sight must be within the lens' Einstein 
radius before significant photometric amplification is observed. Figure |l| gives the image positions 
and amplitudes as a function of u. Note also that the photometric amplification of a quasi-point 
source by a dark lens is achromatic - a fact that is exploited in photometric microlensing searches. 



Center-of-Light Even though the two microlensing images are (by definition) unresolved, they 
are spatially distinct from the nominal source position (Eq. |]), and have non-trivial intensities 
(Eq. ||). We therefore consider astrometry of the center-of-light position in the instance of 
microlensing, which can be performed at several orders of magnitude below imaging resolution 
limits (for examples see Monet et al 1992] , Benedict et al 1995| , |Gatewood &: De Jonge 1995 , Pravdo 



fc Shaklan 1996| , Perryman et al 1997, and references therein). As can be seen from paczyhski 



1996a| Figure 2, the center-of-light is clearly on the symmetry axis between the source and lens. 



Additionally, the lens is in general luminous, and probably unresolved from the images of the 
background source. If we parameterize the relative lens brightness L\ (in general a function of 
wavelength) in units of the unlensed source brightness, relative to the nominal source position, the 
center-of-light is located at: 



Aft, 



center 




(u 2 + 2)(l-L A V^fi) 



u 2 + 2 + L x uVu 2 + 4 
u{u 2 + 3)Vu 2 + 4 



u 2 + 2 u 2 + 2 



L 



x 



where we have written the expression to identify the multiplicative correction (and approximation 
expanded to leading order) to the dark- lens results as a function of L\; in the limit that the lens is 
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dark the expressions simplify to the quantities external to the braces. This center-of-light position 
for a dark lens is included in Figure |l] as a function of u. The observable astrometric displacement 
of the center-of-light from the nominal source position is ( |H0g et al 1995 , Miyamoto Yosh]| 



1995, Paczyhski 1996b, Paczyhski 1995 all give the dark- lens form): 
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(6) 



where we have introduced the angular Einstein radius rE = Re I Di a , and again explicitly 
indicated the dark- lens result and correction (and approximation to leading order in L\) due to 
a luminous lens. In what follows we will assume a dark lens, but will briefly revisit the instance 
of a luminous lens in §Q. Assuming a dark lens, in the microlensing situation described above 
A9 ~ O(rg) is on the order of a milliarcsecond - a value that is within current and envisioned 
astrometric capabilities of ground-based optical interferometers (Palomar Testbed Interferometer 
- PTI (Colavita et al 1994), Keck Interferometer ( |Keck Interferometer Project 1997 ), Very Large 
Telescope Interferometer - VLTI Qvon der Liithe et al 1994| )), space-based global astrometry 
(Space Interferometry Mission - SIM ( Unwin et al 199~7| ), Global Astrometric Interferometer for 
Astrophysics - GAIA ( Lindegren fc Perryman 1996 )), and even filled-aperture ground-based 



CCD differential astrometry flPravdo Sz Shaklan 1996[ ). It is noteworthy that this astrometric 
perturbation signature is at a maximum for u = y/2, and has a value of 2 _3 / 2 re 0.35 te- Note 
also that this perturbation is positive (semi-)definite - the apparent center-of-light is always 
displaced away from the lens. Finally, as in the case of the photometric amplification of a 
quasi-point source, the astrometric displacement is achromatic in the limit of a dark lens, but in 
general is a function of wavelength if the lens is luminous. 



2.2. Microlensing Encounter — Barycentric Observer 

The geometry of a typical, simple microlensing encounter observed in a barycentric (inertial) 
system is depicted in Figure ^. The source, lens, and observer are taken as free to move linearly 
on the time scale of the lensing event. Consequently, viewed in a plane containing the lens and 
normal to the unperturbed source line-of-sight, the relative source-lens trajectory is approximately 
linear during the encounter. Without loss of generality we may define a coordinate system in 
which the source appears stationary and assume that only the lens is in motion. Defining the x-y 
plane of that coordinate system in the transverse lens plane centered on the source (projection), 
and taking the x-coordinate along the trajectory of the lens (see Figure §), we can write the lens 
motion in this system as: 

V\t trnax) 1 
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with v as the relative lens-source transverse speed, b m i n as the minimum source-lens impact 
parameter, and t max the time of (transverse) lens-source closest approach or maximum 
amplification (b s = b m i n at t = t max ). This relative trajectory yields a simple expression for the 
dimensionless source-lens transverse separation u(t) (and its vectorial counterpart u(i)): 




Jb 2 mm + v*(t-t max y 
u{t) = W E = 

u(i) = XL lens (t)/R E = |^| ( ? ) 
where we have defined dimensionless niinimum impact paxeinietGr p = b m i n 

/Re, the characteristic 

time or time scale for the microlensing event to = Re/v - the time required for the lens 
to move transversely one Einstein radius, and a normalized (dimensionless) time coordinate 
[t] = (t — t max )/tQ. To predict the time-dependent microlensing observables, u(t) or u(t) (Eq. 0) 
is inserted into the expressions for the combined image light amplification (Eq. ||) and apparent 
astrometric perturbation (Eq. ^). [Paczynski 1996a| Figure 5 gives sample photometric signature 
lightcurves for simple microlensing events for several values of the minimum impact parameter 
p. This lightcurve shape is precisely the photometric signature for simple microlensing; current 
microlensing survey projects (OGLE - paczynski et al 1995| , MACHO - Alcock et al 1996, EROS 



[Renault et al 1997 , DUO - [Alard et al 1995b| ) first started observing such lightcurves in 1993, 



and by this time have observed many such lightcurves that agree well with this simple theoretical 
expectation. 

As seen by the barycentric observer the source apparently performs an elliptical excursion 
from its unperturbed position QWalker 1995|) . The apparent source astrometric excursion is 
straightforwardly obtained by inserting (the vectorial components of) u(t) (Eq. ^) into A6(u) 
(Eq. §: 

- _ R E -u([t]) _ te j -[t] \ 

and depicted in Figure || for sample values of p (p = 0.2, 0.4, 0.6, 0.8, and 1). It is likewise 
straightforward to derive the eccentricity of this excursion in terms of p: 



p 2 + 2 



The excursions become more eccentric with decreasing p, as the maximum intensity of the two 
images is approximately equal (Eq. The maximum scalar magnitude of the astrometric 
signature, 2 _3//2 te, is the same for all photometrically observable (p < 1) events - it is in fact the 
same for all events with p < y/2. The angular Einstein radius is given in Figure ^| as a function of 
lens distance for bulge and LMC microlensing events and a representative range of lens masses (m 
= 0.05, 0.1, 0.25, 0.5, 1 Mq). For example, a m = 0.1 Mq object at 8 kpc lensing a LMC source 
has an angular Einstein radius of 3 x 10 -4 arcseconds or 300 microarcseconds (1 microarcsecond, 
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Fig. 1. — Image positions and intensities for a dark lens. Left: The two image apparent positions 
and the center-of-light apparent position (in the transverse lens plane) as a function of u. Right: 
image (image 1, image 2, combined) magnifications (dimensionless) as a function of u. 




Fig. 2. — Microlensing encounter geometry. The geometry and differential frame coordinate system 
of a microlensing event is depicted in a plane normal to the source line-of-sight containing the lens. 
The relative source- lens trajectory moves from left to right (direction of increasing x), and is 
parameterized by the lens speed and minimum impact parameter relative to the source. 
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Fig. 3. — Simple microlensing encounter center-of-light motion. Shown here for a range of sample 
p-values are the motions of the apparent center-of-light position traces relative to the nominal 
source position for encounter times -20 < [t] < 20 (Eq. |8|). The profiles are elliptical (note the 
anisotropy of the scale), and are seen to become more eccentric with decreasing p- value (Eq. ||). 
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abbreviated //as, equals 10 6 arcseconds) - and a maximum astrometric signature of roughly 100 
//as. 

It is noteworthy that the time evolution of the astrometric excursion is non-uniform. Far from 
lens-source closest approach the apparent source motion is slow, while near the closest approach 
the source motion is significantly higher. This behavior is depicted in Figure ^ (where points are 
plotted along the excursion trajectory at equal time intervals), and can be seen by differentiating 
the astrometric excursion with respect to time: 



dA6([t]) = -T E j p 2 - [t] 2 + 2 
d [t] (P 2 + W + 2) 2 j 2p[t] 



(10) 



Clearly for [t] 2 S> p 2 (i.e. far from maximum amplification) the excursion rate magnitude goes as 
[t]~ 2 . This behavior is significant in that the astrometric excursion described in Eqs. || and 10 



appears very different than the time-harmonic astrometric excursion expected if the source has 
massive gravitational companions QLestrade et al 1994] , Benedict et al 1995| , patewood 1996| ). 



Finally, it is illustrative to compare the time scales for the photometric and astrometric 
perturbations caused by a microlensing encounter. Figure |5| plots the photometric amplification 
and astrometric perturbation magnitude for positive (normalized) encounter times for a p = 0.4 
microlensing event. The photometric amplification decays from maximum to nominal level in one 
single time scale (to)- In sharp contrast, the astrometric effects are much more persistent, roughly 
requiring factors of 30 and 300 more time to decay to 10% and 1% of the maximum astrometric 
signature respectively. The MACHO collaboration reports a typical time scale to f° r a microlensing 
event to be on the order of one month. Figure || makes the point that measurable microlensing 
astrometric perturbations in events with such a to will span a period of years - depending on the 
astrometric sensitivity. 

The case for astrometric observation of microlensing events is both clear and compelling. 
Barycentric photometric measurements alone constrain the normalized impact parameter p and 
time scale to for a microlensing event through Eq. || flPaczyriski 1996a) Figure 5). Astrometric 
measurements (in lieu of or in addition to photometric measurements) made by a barycentric 
observer additionally constrain the angular Einstein radius Te and the lens transverse motion 
direction (orientation of our x-axis) for the microlensing event through Eq. ^. These two quantities 
taken together are sufficient to compute the (relative) proper motion of the lens object. However, 
because the distance to the lens is not directly established, no direct inference can be made about 
Re, and thereby the mass and transverse velocity of the lens. 



2.3. Microlensing Encounter — Near-Earth Observer 



If the microlensing astrometry is observed by an (near) Earth-based instrument over an 
extended period of time, parallactic effects due to the finite source and lens distances become 
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important as the assumption of linear relative motion in §[T^ is broken by the motion of the 
observer. Figure |6| depicts an example of astrometric trajectories as viewed by barycentric and 
terrestrial observers for an arbitrary p = 0.4 microlensing event. An LMC source is assumed to 
be lensed by a 0.1 Mq lens at 8 kpc (r# ~ 300^/as). It is clear that astrometry sufficient to 
measure the microlensing astrometric excursion will also observe the relative parallactic motions 
of the source and lens. So unlike the case of barycentric observations, the distance to the lens 
(actually, the relative source/lens distance) is accessible by this parallax measurement, and a 
model-independent estimate of the lens mass and (transverse relative) velocity can be derived. An 
independent model for source distance and proper motion allows the lens distance and motion to 
be estimated. A clear case of making virtue from necessity. 

An observer at barycentric 3-position b (measured in AU) observes an object in barycentric 
direction n^, to have an apparent parallactic deflection: 

7r [n 6 (n 6 • b) - b] 

with 7r as the parallax of the source (given in arcseconds by the reciprocal of the distance to the 
source in parsecs). Both the source and lens are at finite distance, so both acquire parallactic 
displacements. In the two-dimensional microlensing coordinate system of § |2.2| , we can write the 
time-dependent correction to the (normalized) source-lens separation u(i) (Eq. [7|): 

Au(t) = [n lens (h lens ■ b(t)) - h(t)} - 7Tsource [h source (h source ■ b(t)) - b(t)] 

rE r E 

« nens ~ 7Wce [n(n ■ b(t)) - b(t)l = — [n(n • b(t)) - b(t)} (11) 

rE r E 

assuming n/ ens « n source n (bnjn ~ O(10~ 9 )), defining the relative lens-source parallax 

n 

= Ttiens ~ 71 source: and with n and b rotated into the microlensing coordinate system. With this 
parallactic correction to u(i), Eq. || predicts the astrometric excursion observed by the terrestrial 
observer (an example of which is depicted in Figure ^). However, it is important to remember 
that Eq. [8| refers to the excursion relative to the unperturbed source position, which itself now 
appears to move with time in an inertial frame due to parallactic effects. This formulation is 
useful, because, as we will argue below, this observational problem lends itself to narrow-angle 
differential astrometry techniques. 

There are several remarkable points concerning the finite distance correction to the 
microlensing astrometry. First, as noted above, fitting a model based on Eq. |8| to such astrometric 
data allows the direct estimate of all the physical parameters for the lens. Strictly speaking, this 
statement assumes D so and the source proper motion can be separately established (or inferred) 
by other means. For instance, the distance to the lens (in parsecs) is simply given by: 

Dlo = Ds ° (12) 

i + d so ii y 1 

with D so in parsecs and II in arcseconds. However, the lens mass is a special case in that it can 
be estimated just from quantities we observe directly, namely rE and n. This fact can be seen by 
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inverting Eq. to solve for m: 



m 



R% 



AG D so x (1 



AG Ds ° rE 1 



x 



— — be] 

ag n [m 



(13) 



where the last equality comes from the fact that x/(l — x) = 1/ D so Il with D so in pc and II in 
arcseconds. The units of te are in radians and II in arcseconds, and the [pc] factor is the length 
of a parsec expressed in the length units of c and G (e.g., meters). The dimensions of Eq. [l3| are 
lacking in elegance, but the result is that we have exchanged the uncertainty in an independent 
estimate of D so for the uncertainty in the definition of an astronomical unit and the quantities that 
can be directly measured in a suitable differential astrometric frame. The second interesting point 
is to note the breaking of the time symmetry around t max by the parallactic correction. Figure ^ 
demonstrates this point in a particular instance as the barycentric t m ax time points are labeled 
on the two excursion trajectories. The degree of symmetry breaking is dictated by the observer's 
orbit phase and source position in the sky. In general this parallactic symmetry breaking leads 
to asymmetric lightcurves and different times for maximum amplification for the barycentric and 
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), and has been observed ( 


Alcock et al 1995). 



3. Astrometric Observations 

During a microlensing encounter photometric observations alone constrain the impact 
parameter p and to through Eq. ||, while microarcsecond-class astrometric observation can be used 
to uniquely determine the fundamental parameters of the lens (mass, distance, and transverse 
velocity) without appealing to a lens population model. Strictly speaking, computing the lens 
distance and transverse velocity requires the source distance and proper motion to be separately 
established, while the lens mass can be estimated without appealing to the source distance 
provided the relative parallax between lens and source can be established (Eq. |b|). Such a 
measurement can be made on the basis of differential astrometry provided an astrometric frame 
in which the background source is quasi-stationary can be established. Currently LMC, SMC, and 
bulge microlensing events are identified by photometric programs searching rich fields of objects 
at roughly the same distance as a lensed source. A differential astrometric frame formed from 
such objects would then have roughly the same parallactic and proper motions as the lensed 
source, and the apparent source excursion could be measured against this reference to establish 
te arid II. This result is compelling, because while microarcsecond wide-angle astrometry requires 
a space-based platform ( Lindegren fc Perryman 1996 , Unwin et al 1997]) , microarcsecond-class 



differential astrometry over narrow fields is possible from the ground ( |5hao fc Colavita 1992| ). 
These considerations strongly suggest a program to perform differential astrometry on microlensing 
candidate events detected in photometric surveys; in such a program the lens mass would be 
directly determined, and the lens distance would be infered based on a model distance to the 
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source. 

It is interesting to investigate the sensitivity that such an astrometric program would 
produce. To address this question we have constructed a simulation code that creates synthetic 
photometric and astrometric measurement sets, and then fits a parametric microlensing model to 
these measurement sets. The measurement synthesis model uses a parametric microlensing model 
(i.e. lens motion relative position angle, p, te, to, t max , nominal source amplitude, II, nominal 
source position, and source proper motion in the frame) and specifications of the measurement 
sequences to simulate (time intervals and measurement frequencies for photometric and astrometric 
sampling, measurement error model - e.g. the sigmas for the zero-mean gaussian errors applied 
to the simulated measurements). Then a microlensing model fitting procedure is used to 
simultaneously fit the synthetic astrometric and photometric datasets and extract estimates of the 
microlensing event parameters. We have studied parameter estimation performance using both 
derivative-based (Marquardt least-squares) and derivative-independent (downhill simplex) fitting 
methods - the results are in good agreement with each other. The simulation code is structured to 
perform this operation in a Monte Carlo mode, so error distributions in the extracted microlensing 
parameters may be empirically derived as a function of physical parameters and measurement 
error models. 

Figure shows sample outputs of the synthetic observation and fitting process for an LMC 
microlensing event. For this example the lens is again assumed to be at distance of 8 kpc (II = 100 
//as), te = 300 //as (Re ~ 2.4 AU, m ~ O.IMq), to = 0.1 yr (relative lens-source transverse speed 
of 120 km s" 1 ), and a lens motion position angle of 30 deg. Operationally we assume the event is 
detected photometrically, and an astrometric campaign starts - this implies astrometry near and 
after maximum amplification. Assuming 3% photometry and 10 //as differential astrometry (the 
sigmas of the zero-mean gaussian errors added to the measurement sequences) referenced to the 
unlensed source position, we fit a parametric microlensing model to the combined data set. The 
fit is seen to reproduce the measurement sets faithfully, and predicts the microlensing parameters 
accurately. Figure || shows the microlensing parameter residuals obtained in 500 instances of the 
event depicted in Figure 0. In particular, we observe 2% and 16% fractional sample standard 
deviations on the r# and II residual distributions respectively. (Fitting Gaussian profiles to the 
central parts of the residual distributions result in error estimates that are ~ 20% better that the 
sample standard deviations). Using the sample standard deviation figures, the error in estimating 
the lens mass is 16% in Eq. [13], clearly dominated by the parallax error. Further, if we were to 
assume a 10% error in estimating the source distance at the LMC, Eq. |l2| yields a lens distance 
accurate to 13%. Such an estimate of the lens distance would unequivocally identify the lens as 
a member of the galactic halo, and exclude the possibility that the lens was either in the galactic 
disk or the LMC itself at high confidence flSahu 1994 Pates et all99^ ). 



By Eq. 13, a given mass lens lies along a particular contour in the space of II vs. r#. Armed 
with our event simulation code we have surveyed this II vs. te phase space of microlensing 
events for experiment sensitivity to microlensing parameters - with a particular emphasis on the 
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lens mass. Table | gives a summary of the set of microlensing and astrometry parameters we 
considered in our Monte Carlo runs. We ran 500 instances of each parameter permutation given 
in the table - a total of 384 cases in all. In addition to the parameters specified in Table [l], 
in each of these experiments we assumed conditions similar to those described above in Figure 
|?]: namely 3% photometry for event detection, and a 5-yr astrometric sequence starting at 
tmax — to j and sampling with uniform period 0.1 to- Microlensing parameter fits are made to the 
astrometric dataset combined with a 3% photometric dataset that spans the interval t max ± 1 yr. 
We considered fitting only the astrometry sequence in a few selected cases, but invariably found 



significantly degraded parameter estimates (similar remarks can be found in H0g et al 1995 ) 



Figure || shows the variation of fractional lens mass error vs. astrometric error for three 
illustrative cases: r# = 100 //as, II = 50 //as (m = 0.02 Mq); r# = 300 //as, II = 100 //as (m 
= 0.1 Mq); and te = 1000 //as, II = 200 //as (m = 0.6 Mq). In each of these cases p = 0.4 
and to = 0.1 yr. Here the lens mass error is estimated from the observed uncertainties (residual 
sample standard deviations) in te and II (and the generally nonzero covariance). This behavior is 
suggestive that the mass error scales by a power law of the astrometric error - we find these cases 
to be typical of the range of Monte Carlo cases considered. 

Ignoring the correlation term, we can estimate the fractional mass error as a function of the 



fractional errors in te and II from Eq. |13 



m 



<4ffL on 

n 2 




P L + H^-m/- + ^ (14) 



where the second equality comes from crudely estimating the uncertainties in and II by a as t- 



Figure 10 shows a scatter plot of the observed fractional mass error from all 384 of our Monte 



Carlo cases against the estimated fractional error given by Eq. 14. The agreement between the 



observed and estimated quantities is good, but we find Eq. |14] overestimates the fractional error 
at large values of a as t ■ A power-law fit to the data indicates that the observed error scales as the 
estimate in Eq. [L4| (hence a a st) to the 0.9 power. We attribute this modest sub-linear scaling to 
the fact that astrometric fits supported by associated photometry resolve te and II slightly better 
than the naive estimate of a as t- We attribute the observed scatter in the data to the correlation 
terms which are included in the observed error estimates, but are neglected in Eq. |14|. We also 
note that the to = 0.2 yr cases on average have lower observed error and smaller scatter than the 
corresponding to = 0.1 yr cases. This is not particularly surprising, as the slower time evolution 
allows for extended observation of the relative parallactic effects. 

Finally, we have argued for and simulated differential astrometric experiments where we 
assume astrometry in a frame where the source is stationary. This reference frame must be 
established by observing field objects near the nominal source position. For LMC, SMC, and bulge 
events where there are many objects at small 5D/D with the source, the reference frame will 
absorb common parallactic motions. However, there will be residual frame drifts and rotations 
resulting from unknown velocity dispersion among the reference objects (and source). The residual 
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Fig. 4. — Microlensing Encounter Angular Einstein Radius. Left: the angular Einstein radius (te) 
for bulge sources as a function of lens distance for lens masses between 0.05 and 1 Mq (m = 
0.05, 0.1, 0.25, 0.5, 1 Mq). Right: the same quantity for LMC events. The astrometric signature 
magnitude for microlensing events is given by 2~ 3 / 2 rg; ~ 0.35 te- Thus a m = 0.1 Mq object at 8 
kpc lensing a LMC source has an angular Einstein radius of 300 //as - and a maximum astrometric 
signature of roughly 100 /ias. 
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Table 1: Microlensing Monte Carlo Parameter Space. This table gives the set of microlensing and 
astrometric accuracy parameters we considered in our Monte Carlo runs. We ran a complete set of 
cases spanning all possible permutations of these parameter values (384 cases in all). 
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Fig. 5. — Comparison of Microlensing Astrometric and Photometric Effects. This figure shows 
the amplification factor and magnitude of the astrometric perturbation as a function of encounter 
time for an event with p = 0.4. The photometric amplification is seen to decay from maximum 
amplification to a nominal value of 1 (no amplification) in roughly 1 event time scale (to). 
By comparison, the astrometric perturbation (normalized to the perturbation at t = t max , 
r E p/(p 2 + 2)) increases to its maximum value at [t] ~ 1, and then decays to zero in a time scale 
several orders of magnitude larger than the photometric amplification. If to ~ 1 month as reported 
for LMC microlensing events by the MACHO collaboration, measurable astrometric deflections in 
such events would last over periods of years - depending on astrometric sensitivity. 
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Fig. 6. — Parallactic Perturbation to Bary centric Microlensing Astrometry. This figure shows the 
astrometric deflection for a p = 0.4 microlensing event as viewed by barycentric and terrestrial 
observers over a period of ± 50 to = ± 5 yr. Marks on the two trajectories are given at to 
= 0.1 yr intervals. (The parallactic motion of the background source - 20 fias - is removed to 
facilitate comparison.) The background source is taken in the LMC (D so ~ 50 kpc), and the lens 
is taken to be in the galactic halo (D[ ~ 8 kpc), with an angular Einstein radius of 300 ^as (m ~ 
0.1 Mq). Astrometry sufficient to measure the microlensing perturbation will also measure the 
parallax effects. A parallax measurement can be used to estimate the relative source/lens distance 
to the lens, which in turn allows the lens mass, distance, and relative transverse velocity to be 
inferred for individual events. The time marks corresponding to (barycentric) t max are shown. 
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Fig. 7. — Sample Microlensing Model Fitting. Here we show an example of fitting a microlensing 
model to synthetic terrestrial photometry and astrometry datasets for a microlensing encounter. 
The critical parameters for the event are a lens motion position angle of 30 deg, p = 0.4, te = 
300 //as, and II = 100 /ias (m = 0.1 Mq). We assume the event is identified photometrically, 
and differential astrometric measurements commence after that detection. The microlensing model 
described in §[2] was simultaneously fit to both the photometric and astrometric data. Shown in 
each are the simulated data, true values, and the model fit. Left: the photometric lightcurve 
results. In this example we assume 3% RMS error photometry. The time units on the x-axis 
are plotted relative to the barycentric t max . Once a microlensing interpretation seems likely, we 
commence astrometric measurement - this time point is shown. Note that the time of maximum 
amplification for the terrestrial observer is offset relative to a barycentric observer. Further, the 
lightcurve is slightly asymmetric with respect to the time of maximum amplification. Both of 
these effects are well known for terrestrial microlensing observation. Right: the corresponding 
depiction for the astrometry sequence relative to the nominal source position. The simulated 
10 //as RMS error astrometric measurements begin shortly before maximum magnification, and 
continue for 50 to after maximum magnification. The true excursion trajectory is shown over the 
complete excursion, but the microlensing model fit prediction is rendered only for the interval of 
the astrometric measurements. 
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Fig. 8. — Sample Monte Carlo Fit Residual Distributions. Here we show a sample of the 
microlensing model fit residual distributions obtained in 500 measurement sequence instances 
taken from the microlensing event depicted in Figure f^. Measurement parameters are the same 
as discussed in Figure [?]: 3% photometry and 10 /uas differential astrometry. Distributions are 
shown for lens motion position angle (top left), normalized impact parameter (top right), angular 
Einstein radius (lower left), and relative parallax (lower right). These astrometric experiments are 
seen to yield fractional uncertainties in the angular Einstein radius and parallax of 2% and 16%, 



respectively, yielding a 16% error in the lens mass estimate by Eq. 13. Further, combined with 
a separate model for the source distance (assumed 10% accuracy), the relative parallax estimate 
preformance yields a 13% estimate for the lens distance, easily sufficient to unequivocally establish 
its location in the halo. 
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Fig. 9. — Fractional Mass Error vs. Astrometric Error: Three Examples. Here we show three 
samples of fractional mass error vs. astrometric error from our ensemble of Monte Carlo results. 
The three examples are: = 100 //as, II = 50 //as (m = 0.02 Mq); te = 300 //as, II = 100 //as 
(m = 0.1 Mq); and rE = 1000 /tas, II = 200 //as (m = 0.6 Mq). Not surprisingly, we find the 
fractional mass error to scale as a power law of the astrometric error in all our Monte Carlo results. 
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Fig. 10. — Scatter of Observed Mass Error vs. Estimated Mass Error. Here we show the scatter 
of the observed fractional mass error vs. the estimated fractional mass error given in Eq. 14. Each 
data point represents the fractional mass error obtained in 500 instances of measurement sequence 
fitting. Data sets for to = 0.1 yr, 0.2 yr and p = 0.4, 0.8 are shown separately - the intrinsic 
scatter in the to = 0.2 yr results is noticeably smaller. The data are seen to correlate reasonably 
well with the crude mass estimate of Eq. [ij, but it overestimates the mass error at larger values. 
A simple parametric fit shows the observed fractional mass error to scale (with moderate scatter) 
to the estimated fractional mass error to the 0.9 power. We attribute the scatter in the data to the 
correlation in errors included in the observed error calculation but explicitly ignored in Eq. 14. 
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linear frame drifts are large (mas yr _1 ) on the scale of the microlensing astrometric excursion. 
The frame rotations will be small (O(10 -4 ) rad yr _1 ), but to the extent that the source is not at 
the center of the astrometric frame this too adds an effective linear drift to the differential frame 
of order 100 //as yr _1 . So in practice the microlensing fit process must account for linear frame 
drifts that are large compared with the scale of the microlensing excursion. It is straightforward 
to extend the microlensing parameter model to include a frame drift. Figure 11 gives an example 
of allowing for a random drift in the astrometric frame, and using an extended microlensing model 
to solve for this quantity. The physical parameters of the microlensing event in Figure [ll| are the 
same as those used in Figure 0: te = 300 //as, II = 100 //as, to = 0.1 yr; and we have taken 50 
//as astrometry. To this we have added 1.5 mas yr -1 of frame drift in a random orientation to the 
lens motion. The microlensing fit does an acceptable job of identifying the frame drift, bolstered 
by the time base of astrometry at late encounter times (not shown in the figure). We do not find 
a coupling between frame drift and the quantities of interest in determining the lens mass: r# and 
II. We have added random frame drifts to several of the Monte Carlo cases described in Table [E], 
and find that such drifts do not significantly effect our simulation results. This is reasonable, as 
both te and II are estimated from the curvature of the astrometric path. 



4. Summary and Discussion 

We have specialized the Miyamoto and Yoshi suggestion to measure the astrometry of the 
lensing event photocenter. We agree with their conclusions: that high-precision astrometric 
observation of MACHO microlensing events hold the promise of determining fundamental lens 
parameters (mass, proper motion, transverse velocity) in a model-independent fashion. We 
maintain that it is sufficient to measure the motion of the center-of-light. In particular, in 
many cases accurate differential astrometry is sufficient obtain the lens mass without additional 



assumptions (in slight contrast to the earlier suggestion by Walker 1995 , who ignored the value 
of parallactic effects), and reasonable lens distance and transverse velocity estimates can be 
obtained from an independent model of source distance and proper motion. Such narrow-angle 



differential astrometry is possible from the ground (Shao & Colavita 1992). Alternatively, 
wide-angle microarcsecond-class astrometry can simultaneously determine the source and lens 
position and kinematic parameters without external information, and again the lens mass. Clearly 
the potential for probing the physical parameters of the putative halo object population by 
astrometric techniques is enormous. 

A program to probe microlensing events photometrically detected in the galactic bulge seems 



plausible for the planned Keck Interferometer (KI - Keck Interferometer Project 1997 ). KI 
requires a bright guide star to track atmospheric fluctuations of the interferometric fringes. The 
brightest bulge objects are 16th magnitude, within the fringe tracking capabilities for the two 10 m 
apertures. The expected 10 - 20 //as astrometric performance of KI yields microlensing parameter 
estimates sufficient to constrain lens mass and distance parameters for individual events, which 
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Fig. 11. — Microlensing Astrometric Fit With Frame Drift. Shown is the result of an astrometric 
experiment where we have added a linear drift into the differential astrometric frame. The 
microlensing parameters of this event are the same as used in Figure [?]: ve = 300 //as, II = 
100 //as, to = 0.1 yr; and we have taken 50 //as astrometry. We added 1.5 mas yr _1 of frame drift 
in a random orientation to the lens motion (vectors indicating the frame drift and lens motion 
direction are shown). The microlensing fit does an acceptable job of identifying the frame drift. 
We do not find a coupling between frame drift and the quantities of interest in determining the lens 
mass: te and II. We have added random frame drifts to several of the Monte Carlo cases described 
in Table |l|, and find that such drifts do not effect our simulation results. 
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will give profound insight into the nature of these objects. 

However, events in the LMC or SMC are not detectable from the Keck site, both because 
of geography and sensitivity. The brightest objects in the LMC are 17-18 magnitude, arguably 
fainter than the tracking capabilities of the KI. However, the declination of these fields would 
require KI zenith angles that severely degrade the astrometric performance. A large aperture 
astrometric interferometer in the Southern hemisphere such as the VLT interferometer (VLTI - 



von der Liithe et al 1994 ) could measure Magellenic cloud events, and in particular determine the 
lens mass and distance with sufficient accuracy to resolve many of the current issues regarding 
their nature. Such measurements would be very challenging but are clearly very compelling. 

A number of authors have specifically suggested the application of planned space-based global 



astrometric techniques to analyze these events (H0g et al 1995, Miyamoto fc Yoshi 1995, Paczyhski 



1998). In space-based applications astrometric references can be drawn from a global astrometric 
frame tied to extragalactic objects, and the positions, proper motions, and parallaxes of these 
references are known to a few microarcseconds in a quasi-static frame. Thus the necessity 
of establishing a narrow-angle relative frame for differential astrometry is removed. Further, 
astrometry at late times identifies the source proper motion and parallax in the global frame - 
thereby establishing the source motion and distance. With the source distance and kinematics 
established, the lens parameters are all uniquely determined. Assuming ground-based observations 
measure lens masses and distance in the manner we describe, we believe the role of space-based 
microlensing astrometry programs by planned astrometric space missions such as SIM (|U nwin etj 



al 1997 ) and GAIA ( [Lindegren fc Perryman 1996 ) will be in probing the precision positions and 



particularly kinematics of the lensing objects. If ground-based measurements of Magelenic cloud 
events are not possible, then SIM and GAIA seem well-suited to offer definitive answers on the 
nature of the lenses. 

Finally, in the near future it is possible that CCD-based astrometry could make detections of 
microlensing astrometric perturbations, and possibly make rough estimates of lensing parameters, 
and/or breaking some of the degeneracies in photometric microlensing observations (see below). 
Pravdo and Shaklan ( Pravdo &: Shaklan 1996] ) report night-to-night astrometric repeatability of 



200 fias in data taken at the Palomar 5 m telescope, and speculate that limits might approach 
100 /ias at the 10 m Keck Telescope. In further assessing these prospects we anxiously await the 
results of several nights of Keck observations recently made by Pravdo and Shaklan ( Shaklaij 



19971) . 



One of the key assumptions we have made in this work is the assumption of a dark lens. 
This assumption is plausible given the success photometric programs have had in fitting dark lens 
amplification models to photometric data. However, the instance of a luminous lens is possible and 
interesting ( Miralda-Escude 1996| , paczyhski 1996b ), and the astrometric model derived here can 



be augmented in a straightforward way. Figure [T^ shows an example of a fit to a dataset generated 
with a luminous lens model unresolved from the lensed source. The physical parameters (bulge 
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event, 10 kps source / 5 kpc lens distance, 0.1 Mq lens, lens/source brightness ratio of 0.1) were 
selected to be comparable to the example given in Figure [?]. The model fit faithfully reconstructs 
the input physical parameters, including the correct attribution of source and lens brightness. 
One operational question that arises is how does one determine whether a dark lens model is 
appropriate for a given microlensing event. While one could straightforwardly test the luminous 
lens hypothesis by adding a relative source/lens luminosity parameter to the microlensing fit 
model as presented here, a more obvious and compelling resolution to this question is contained 
in the possible chromaticity of the astrometric observables as described in Eq. ^. If the lens is 
luminous, then its spectral content is in general different from that of the source - implying that 
both the photometric and astrometric observables will be functions of wavelength. Making the 
observations in a variety of spectral bands will identify the relative source- lens intensity and color, 
and provide the necessary data to robustly extend the microlensing model to the general case of 
luminous lenses. We defer a more systematic analysis of the luminous lens case to future work. 

The closely related problem of image blending is discussed in recent work by Alard ( Alard| 
1996), and Wozniak & Paczyhski ( Wozniak fc Paczyhski 1997] ), who consider the possibility 



that a second source (possibly unrelated to either source or lens) is unresolved from the image. 
Wozniak & Paczyhski find that the degeneracies in photometric observation of such events result 
in systematic errors in estimating lensing parameters. Based on our preliminary successes in 
correctly distinguishing source and lens luminosities, we concur with the speculation put forward 
by Wozniak & Paczyhski that multi-spectral astrometry and photometry breaks the degeneracy 
in (some subset of) blended events, and point this out as a particularly important case for future 
study. 



Astrometric Detections — Non-MACHO Events A number of authors have suggested to 
broaden the applicability of the astrometric techniques to generic microlensing events ( [Hosokawa 



et al 1993| , |Miralda-Escude 1996 , Paczyhski 1996b, Paczyhski 1998). These events could potentially 
be detected astrometrically in programs that concentrated on high proper motion objects (so as to 
sweep-out larger solid angles), or as a part of broader companion search program (something we 
have integrated into our PTI program - |Colavita et al 1994 ). While much of the phenomenology 
we have developed in §|2] is directly applicable, there is a practical difficulty in establishing lensing 
parameters in such events by differential means. The first is the absence of a ready supply of 
reference objects that share common parallactic motions as the source. The fact that the rich 
LMC, SMC, and bulge fields used in the photometric surveys naturally yield an abundance of 
reference objects for which 5D/D is small makes events in these fields unique. Without such a 
common parallactic reference, the systematic errors in the determination of II will be too large to 
establish a precise lens mass from ground-based differential astrometry. Such events would seem 
to be best studied by space-based, global astrometric techniques. 
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Fig. 12. — Luminous Lens Microlensing Model Fitting. Here we show an instance of fitting a 
microlensing model with an unresolved luminous lens to synthetic terrestrial photometry and 
astrometry datasets for a microlensing encounter similar to that shown in Figure 0. The parameters 
for the event are a lens motion position angle of 30 deg, p = 0.4, = 290 //as, and IT. = 100 //as 
(m = 0.1 Mq) and a lens/source brightness ratio of 0.1. Again, we assume the event is identified 
photometrically, and differential astrometric measurements commence after that detection. The 
microlensing model described in §[2] including a parameter for the relative lens/source brightness 
was simultaneously fit to both the photometric and astrometric data. Shown in each are the 
simulated data, true values, and the model fit. Left: the photometric lightcurve results. In this 
example we assume 3% photometry error. The time units on the x-axis are plotted relative to 
the barycentric t max . Right: the corresponding depiction for the astrometry sequence relative to 
the nominal source position. The simulated 20 //as astrometric measurements begin shortly before 
maximum magnification, and continue for 30 to after maximum magnification. The fit is seen to 
faithfully reproduce the simulated datasets, and converge to the input model values including the 
appropriate source and lens brightness, even in the presence of frame drift. 
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Complex Lenses While the majority of photometrically detected events are consistent with 
single lens hypothesis, a number of binary lens events have been reported ( [Udalski et al 1994 , Alard 
|et al 1995a , Bennett et al 1996| ) . In a recent preprint Dominik (Dominik 1997) argues that 
photometry alone does not uniquely constrain the binary lens parameters. We speculate that 
additional astrometric information would break the degeneracies among various hypotheses in 
binary lens events through straightforward extensions of the astrometric theory developed here. 
We defer the analysis of the binary lens case to future work. 



We gratefully acknowledge the many helpful comments made by the anonymous reviewer, and 
Prof. B. Paczyhski. We also thank Prof. D. Peterson for frequent and unfailing "encouragement". 
The work described in this paper was performed at the Jet Propulsion Laboratory, California 
Institute of Technology under contract with the National Aeronautics and Space Administration. 



REFERENCES 

Alard, C, Mao, S., Guibert, J. 1995, A&A 300, 17 ( |astro-ph/9506T0ll ). 
Alard, C. et al (the DUO collaboration) 1995, Msngr 80, 31. 
Alard, C. 1996, A&A in press ( b.stro-ph/9609165|) . 

Alcock, C. et al (the MACHO collaboration) 1995, ApJ 454, L125 (|astro-ph/9506114|) . 
Alcock, C. et al (the MACHO collaboration) 1997, ApJ 486, 697 ( hstro-ph/9606165| ). 
Benedict, G.F. et al 1995, BAAS 187, #70.11. 

Bennett, D.P. et al (the MACHO collaboration) 1996, Nucl. Phys. Proc. Suppl. 51B, 152 
( |as¥o-ph/9606012D . 

Buchalter, A. and Kamionkowski, M. 1996, ApJ 469, 676 ( |astro-ph/9604144j ). 

Colavita, M.M. et al 1994, Proc. SPIE 2200, 89. 



Dominik, M. 1997, preprint (|astro-ph/9703003D . 

Gates, E.I., Gyuk, G., Holder, G.P, and Turner, M.S. 1997, ApJ submitted (|astro-ph/97TTTTo|) . 
Gatewood, G., and De Jonge, J.K. 1995, ApJ 450, 364. 
Gatewood, G. 1996, BAAS 188, #40.11. 

Gaudi, B., and Gould, A. 1997, ApJ 477, 152 ( |astro-ph/960T030| ). 
Gould, A. 1992, ApJ 392, 442. 



-27- 



Gould, A. 1996, PASP 108, 465G ( |astro-ph/96040ll| ). 

H0g, E., Novikov, I.D., and Polnarev, A.G. 1995, A&A 294, 287. 

Hosokawa, M., Ohnishi, K., Fukushima, T., and Takeuti, M. 1993, A&A 278, L27. 



Keck Interferometer Project 1997 http: / /huey.jpl. nasa.gov/keck 



Lestrade, J.F., Jones, D., Preston, R., and Phillips, R. 1994, Ap&SS 212, 251. 
Lindegren, L. and Perryman, M.A.C. 1996, A&AS 116, 579L. 

von der Liithe, O., Ferrand, D., Koehler, B., Neng-hong, Z., and Reinheimer, T. 1994, 
Proc. SPIE 2200, 168. 

Miyamoto, M. and Yoshi, Y. 1995, AJ 110, 1427. 

Miralda-Escude, J. 1996, ApJ 470, L113 ( |astro-ph/9605138| ). 

Monet, D.B. et al 1992, AJ 103, 638. 

Paczyiiski, B. 1986a, ApJ 301, 503. 

Paczyiiski, B. 1986b, ApJ 304, 1. 

Paczyhski, B. et al (the OGLE collaboration) 1995, BAAS 187, #14.07 
Paczyiiski, B. 1996a, ARA&A 34, 419 festro-ph/9604011| ). 
Paczyhski, B. 1996b, Acta Ast. submitted ( astro-ph /9606060|) . 
Paczyhski, B. 1998, ApJ 494, L23 ( |astro-ph/9708T55|) . 
Perryman, M.A.C. et al 1997, A&A 323. L49. 
Pravdo, S.H. and Shaklan, S.B. 1996, ApJ 465, 264. 
Refsdal, S. 1964, MNRAS 128, 295. 

Renault, C. et al (the EROS collaboration) 1997, A&A 324, L69 ( |astro-ph/9612102i ). 
Sahu, K.C. 1994, PASP 106, 942 ( |astro-ph/9408047| ). 
Shaklan, S.B. 1997, private communication. 
Shao, M., and Colavita, M. 1992, A&A 262, 353. 



Udalski, A. et al (the OGLE collaboration) 1994, ApJ 436, L103 (|astro-ph/9407084|) . 
Unwin, S., Boden, A., and Shao, M. 1997, Proc. STAIF, AIP Conf. Proc. 387, 63. 



-28- 



Walker, M.A. 1995, ApJ 453, 37. 

Wozniak, P. and Paczyriski, B. 1997, ApJ 487, 55 ( |astro-ph/9702194| ). 
Zhao, H-S. 1997, MNRAS submitted ( |astro-ph/9703097j ). 



This preprint was prepared with the AAS IATgX macros v4.0. 



